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Abstract 

The  optimal  spectral  decomposition  (OSD)  method  is  used  to  reconstruct  seasonal  variability  of  the  Black  Sea  horizontally 
averaged  chlorophyll-a  concentration  from  data  collected  during  the  NATO  SfP-97 1818  Black  Sea  Project  in  1 980-1995.  During 
the  reconstruction,  quality  control  is  conducted  to  reduce  errors  caused  by  measurement  accuracy,  sampling  strategy,  and  irregular 
data  distribution  in  space  and  time.  A  bi-modal  structure  with  winter/spring  (February-March)  and  fall  (September-October) 
blooms  is  uniquely  detected  and  accurately  documented.  The  chlorophyll-a  enriched  zone  rises  to  1 5  m  depth  in  winter  and  June, 
and  deepens  to  40  m  in  April  and  35  m  in  August.  The  June  rise  of  the  chlorophyll-a  enriched  zone  is  accompanying  by  near- 
continuous  reduction  of  upper  layer  maximum  chlorophyll-^  concentration. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Observational  (in-situ  and  remote  sensing)  and 
modeling  studies  show  rapid  change  and  deterioration 
of  the  Black  Sea  ecosystem.  The  abruptly  increased 
nutrients  and  contaminants,  that  are  delivered  to  the 
northwest  shelf  zone  with  river  runoffs  and  redis¬ 
tributed  within  the  entire  water  body,  have  induced  a 
severe  basin-scale  eutrophication  of  the  Black  Sea 
water  (Mee,  1992;  Zaitsev  and  Mamaev,  1997;  Yilmaz 
et  al.,  2001;  Oguz  et  al.,  2002a). 
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The  chlorophyll-a  concentration  is  considered  to 
be  the  main  measure  of  phytoplankton  biomass  and 
has  been  intensively  discussed  as  applied  to  the  Black 
Sea  in  a  variety  of  aspects  including  patterns  of  spatial 
distribution,  seasonal  and  inter-annual  variability 
(Vedernikov  and  Demidov,  1993,  1997;  Oguz  et  al., 
2002b;  Kopelevich  et  al.,  2002;  Yunev  et  al.,  2002; 
Moncheva,  2003). 

The  chlorophyll-a  data  come  from  two  major 
sources:  in-situ  measurements  and  satellite  (CZCS 
and  SeaWiFS)  images.  The  in-situ  data  are  noisy 
(noise-to-signal  ratio  up  to  60%),  sparse  and  irregular 
in  space  and  time  (Vedernikov  and  Demidov,  1997). 
The  satellite  data  (SeaWiFS  images)  are  also  noisy 
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with  the  noise-to-signal  ratio  near  35%  for  open  water 
and  higher  for  turbid  coastal  water  such  as  the  Black 
Sea  (Srokosz,  2000;  Oguz  et  al.,  2002b). 

Different  types  of  seasonal  variability  of  the  basin- 
scale  Black  Sea  chlorophyll-a  concentration  have 
been  identified  using  the  two  kinds  of  data  sources 
and  traditional  statistical  methods.  A  bi-modal  struc¬ 
ture  with  winter/spring  and  fall  blooms  is  considered 
to  exist  in  the  phytoplankton  annual  cycle  in  the 
middle  latitudes  (Longhurst,  1995;  Sathyendranath  et 
al.,  1995)  and  must  appear  in  the  Black  Sea  as  well. 
The  February-March  bloom  presents  a  pronounced 
signal  in  the  chlorophyll-a  measurements,  and  so  it 
has  been  clearly  detected  using  both  in- situ  and 
satellite  data  (Vedernikov  and  Demidov,  1997;  Oguz 
et  al.,  2002b,  2003,  among  others).  It  is  not  the  case 
for  the  fall  bloom,  the  existence  of  which  was 
surmised  from  the  satellite  data  analysis  (Oguz  et 
al.,  2002b,  2003),  but  was  not  evidently  proved 
through  the  in-situ  measurements  (Vedernikov  and 
Demidov,  1997;  Yunev  et  al.,  2002). 


A  question  arises:  How  many  blooms  are  there  in 
the  seasonal  variability  of  the  Black  Sea  chlorophyll-a 
concentration?  Since  the  satellite  data  are  more 
regularly  sampled  in  space  and  time  than  the  in-situ 
data,  a  reasonable  guess  may  favor  the  two  blooms, 
but  ultimate  confirmation  of  the  two-bloom  structure 
still  needs  the  in-situ  data  with  its  vertical  sampling 
(subsurface  information).  Since  the  in-situ  data  are 
highly  noisy  and  irregularly  sampled,  the  traditional 
statistical  methods  might  not  be  sufficient.  Methods 
with  the  capacity  to  process  noisy  and  sparse  data 
should  be  used. 

Recently,  an  optimal  spectral  decomposition 
(OSD)  method  was  developed  to  reconstruct  noisy 
and  sparse  scalar  fields  as  well  as  Eulerian  and 
Lagrangian  velocity  data.  This  method  is  robust  to 
various  noise-to-signal  ratios,  number  of  observa¬ 
tions,  and  sampling  strategy  (Chu,  1999;  Chu  et  al., 
2003  a, b).  For  example,  reliable  data  sets  were 
obtained  for  radionuclide  pollution  of  the  Black, 
White,  and  Kara  Seas  using  the  OSD  method 
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Fig.  1.  The  Black  Sea  geography,  bottom  topography,  and  circulations  (Oguz  et  al.,  1993). 
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(Ivanov  et  al.,  2001,  2004;  Eremeev  et  al.,  1994, 
1995;  Danilov  et  al.,  1996;  Ivanov  and  Margolina, 
1996). 

In  this  study,  the  Black  Sea  chlorophyll-a  concen¬ 
tration  data  collected  during  the  NATO  SfP-971818 
Project  (in-situ  data)  are  analyzed  using  the  OSD 
method.  The  main  purpose  is  to  identify  the  basin- 
scale  seasonal  variability  (one  or  two  blooms?).  The 


rest  of  the  paper  is  outlined  as  follows.  Sections  2  and 
3  describe  the  Black  Sea  oceanography  and  features  of 
the  in-situ  chlorophyll-a  observations.  Section  4 
depicts  the  framework  of  the  OSD  method  and 
discusses  its  merits  and  limitations.  Section  5  dis¬ 
cusses  the  reconstructed  chlorophyll-a  field. 

In  Section  6,  reliability  of  the  reconstruction  results 
is  analyzed.  The  main  conclusions  are  presented  in 


16  18  20  22  24 

Salinity  (2),  °/00 

12  14  16  18  20 

Density  (3),  kg  m'3 


Winter 


-0.05  0  0.05  OA 

Salinity  gradient  (2),  °/00  m'1 


0  0.02  0.04  0.06 

Density  gradient  (3),  kg  m'4 

Summer 


Temperature  gradient  (1),  °C  m'1 


18  20  22  24 

Salinity  (2),  °/00 

10  15  20 

Density  (3),  kg  m'3 


-0.02  0  0.02  0.04  0.06 

Salinity  gradient  (2),  °/00  m'1 
0  0.05  0.1 

Density  gradient  (3),  kg  m'4 


Fig.  2.  The  Black  Sea  temperature,  salinity,  and  potential  density  profiles  and  their  vertical  gradients:  (a)  winter  profiles,  (b)  winter  vertical 
gradients,  (c)  summer  profiles,  and  (d)  summer  vertical  gradients. 
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Section  7.  Appendix  A  provides  discussion  on  the 
capability  of  the  OSD  to  fill  considerable  spatial  gaps 
in  observation  coverage. 

2.  Black  Sea  oceanography 

The  Black  Sea  is  the  largest  inland  basin  in  the 
middle  latitudes  with  coastal  and  deep  parts  pro¬ 
foundly  dissimilar  to  each  other  (Fig.  1).  The  deep 
basin  with  a  relatively  flat  bottom  (depth  between 
2000  and  2300  m)  occupies  more  than  60%  of  the 
total  area.  The  wide  northwestern  continental  shelf 
with  depth  about  50  m  takes  most  of  the  remaining 
area.  Thus,  interactions  between  coastal  and  deep 
ecosystems  are  inherent  for  the  Black  Sea.  The  basin- 
scale  circulation  appears  as  a  cyclonic  large-scale 
flow,  the  Rim  Current.  This  jet-current  is  associated 
with  two  principal  cyclonic  gyres  situated  in  the 
western  and  eastern  parts  of  the  basin,  and  anti- 
cyclonic  eddies  dominating  the  coastal  zone  (Fig.  1). 

Hydrological  balance  of  the  Black  Sea  determines 
its  hydrochemical  structure.  The  Black  Sea  has  a  two- 
layer  hydrographic  structure.  Cold  and  fresh  surface 
layer  water,  which  is  controlled  mainly  by  the  riverine 
input,  overlies  warm  and  salt  deep  water,  which  is 
under  the  influence  of  inflowing  salty  Mediterranean 
water  through  the  Bosphorus  Strait.  Thus,  the  density 
of  subsurface  water  is  determined  mainly  by  the 
salinity.  The  halocline  and  pycnocline  coincide  at  a 
depth  between  100  (near  the  centers  of  the  main 
cyclonic  gyres)  and  200  m  (around  the  coastal 
margins). 

Typical  temperature  and  salinity  profiles  (Blatov  et 
al.,  1984)  are  shown  in  Fig.  2.  The  depth  of  the 
pycnocline  overlaps  with  the  lower  boundary  of  the 
Cold  Intermediate  Layer  extending  over  the  entire 
deep-water  area  of  the  Black  Sea.  This  layer  has  a 
temperature  minimum  (6. 5-7. 5  °C)  between  25  and 
150  m  and  considerably  affects  the  hydrological, 
hydrochemical,  and  biological  processes.  The  salinity 
stratification  substantially  governs  the  vertical 
exchange  and  as  a  result,  causes  the  coincidence  of 
the  oxycline  and  chemocline  at  a  depth  between  100 
and  200  m.  The  permanent  presence  of  intense  density 
stratification  facilitates  the  formation  of  anoxic  water 
masses  below  <70=16.1-16.2  kg  m-3.  On  the  other 
hand,  the  upper  mixed  layer  of  low  salinity  responds 


strongly  to  seasonal  heating  and  cooling  and,  as  a 
result,  the  temperature  predominantly  determines  the 
density  in  the  summer  season  (Fig.  2b).  From  spring  to 
fall,  a  shallow  thermocline  forms  at  a  depth  of  about  30 
m  that  is  the  upper  boundary  of  the  Cold  Intermediate 
Layer.  Obviously,  the  above  description  of  the  Black 
Sea  is  far  from  comprehensive.  For  more  detailed 
information,  the  interested  readers  are  referred  to  a 
more  complete  synthesis  of  the  Black  Sea  ocean¬ 
ography  (Murray,  1991;  Ozsoy  and  Unluata,  1997). 


Table  1 

Sources  of  the  utilized  chlorophyll-a  measurements 


Month 

Year 

RV,  cruise 

Overall 
sampling  size 

Method 

January 

1986 

PV,  20 

18 

FI 

1988 

PV,  25 

46 

FI 

1992 

PV,  36 

5 

FI 

February 

1992 

PV,  36 

3 

FI 

March 

1986 

PV,  20b 

45 

FI 

1988 

ML,  49 

80 

FI 

1995 

PK,  33 

13 

FI 

April 

1988 

PK,  18 

33 

FI 

1989 

PV,  28 

4 

FI 

1993 

PK,  30 

10 

FI 

1993 

PV,  41 

37 

FI 

1995 

PK,  33 

13 

FI 

May 

1982 

PV,  12 

15 

Sp 

1988 

PK,  18 

1 

FI 

June 

1989 

PV,  28 

36 

FI 

1991 

PV,  33 

44 

FI 

July 

1982 

AT,  03/1 

66 

FI 

1985 

PV,  19 

9 

FI 

1992 

PV,  37 

23 

FI 

August 

1980 

PV,  9 

21 

Sp 

1982 

AT,  03/1 

14 

FI 

1990 

PV,  32 

9 

FI 

1992 

PV,  38 

7 

FI 

September 

1980 

PV,  9 

59 

Sp 

1983 

PK,  7 

24 

FI 

1984 

ML,  43-1 

32 

FI 

1985 

PK,  12-2 

81 

Sp 

1990 

ML,  53a 

16 

FI 

1991 

PK,  28 

36 

FI 

October 

1990 

ML,  53a 

65 

FI 

1991 

PK,  28 

10 

FI 

November 

1983 

PK,  7 

43 

FI 

1989 

ML,  51 

21 

FI 

1991 

PV,  35 

25 

FI 

December 

1987 

PV,  25 

16 

FI 

1989 

ML,  51 

14 

FI 

1994 

PK,  32 

71 

FI 
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3.  Chlorophyll-0  concentration  data 

Table  1  lists  the  in-situ  chlorophyll-0  measurements 
(1980-1995)  from  the  NATO  SfP-971818  ODBMS 
Black  Sea  Project.  The  data  were  collected  during  26 
expeditions  of  research  vessels  “Mikhail  Lomonosov” 
(ML),  “Ay-Todor”  (AT),  “Professor  Kolesnikov”  (PK) 
and  “Professor  Vodyanitsky”  (PV)  of  the  Ukrainian 
National  Academy  of  Sciences  (Marine  Hydrophysical 
Institute  and  Institute  of  Biology  of  the  Southern  Seas, 
Sevastopol,  Ukraine).  The  chlorophyll-0  concentration 
was  measured  using  fluorometric  (FI)  and  spectropho- 
tometric  (Sp)  techniques.  The  accuracy  of  these 
measurements  will  be  discussed  in  the  next  section. 

The  data  set  contains  about  1100  chlorophyll-0 
concentration  profiles  from  the  surface  to  100  m  depth. 
They  are  sparse  and  irregular  in  time  and  space  (Figs.  3 


Fig.  4.  Time  (by  month)  and  depth  distribution  of  the  observed 
chlorophyll-^  observations. 
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Fig.  3.  Distribution  of  chlorophyll-a  observational  stations:  (a)  winter,  (b)  spring,  (c)  summer,  and  (d)  fall. 
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and  4).  Certain  periods  and  areas  are  well  sampled 
while  others  lack  enough  observations  to  gain  any 
statistically  meaningful  estimates.  For  example,  Octo¬ 
ber  and  December  have  much  fewer  data  points  than 
the  other  months  (Fig.  4).  Even  the  months  with  a 
relatively  large  number  of  data  have  rather  irregular 
sampling  distribution  in  the  horizontal  and  vertical. 

4.  Methodology 

4.1.  Removal  of  inter-annual  trend 

Since  the  main  goal  here  is  to  estimate  the  seasonal 
variability  of  chlorophyll-a  concentration,  the  inter¬ 
annual  trend  ctrend(0  should  be  removed  prior  to  the 
reconstruction, 

c(x,Z,  t)  =  Cobs(x,Z,  t)  ^trend^)?  (1) 


using  the  generalized  Fourier  series  (Ivanov  et  al., 

2001) 

M 

c(X)ZjC)t')  =  A(){zk,t)  H-  ^  ^JAm(<zkfs)Wm(x1zkf 

m=  1 

xei?(zi)  (3) 

where  M  is  the  truncated  mode  number,  *Fm(x,  zk)  and 
Am(zk,  t)  are  the  orthogonal  basis  functions  (or  called 
modes)  and  the  spectral  coefficients,  respectively; 
R(zk )  is  the  area  bounded  by  the  Black  Sea  lateral 
boundary  P(zk )  at  depth  zk. 

The  eigenfunctions  {F m(x,  zk)}  of  the  horizontal 
Laplace  operator  with  the  basin  geometry  and 
homogeneous  Neumann  boundary  condition  at  T(z) 
are  taken  as  the  basis  functions, 

V^m=-4*Pmi  n-V/,¥gr  =  0, 

m=\,2,...,M,  (4) 


where  cobs(x,  z,  t )  is  the  measured  chlorophyll-a 
concentration  at  location  (x,  z)  and  time  t.  The  Black 
Sea  chlorophyll-a  concentration  does  not  have  any 
evident  interannual  trend  before  1988,  and  has  distinct 
interannual  trends  between  1988-1992  and  1992- 
1995  (Yunev  et  al.,  2002).  The  traditional  approach 
(Saltzman,  1983)  is  used  to  estimate  the  interannual 
trends  for  the  two  periods  1988-1992  and  1992-1995, 

ao  +  ait,  fc[1988,1992]  ,  . 

bo  +  bif  te[1992, 1995]’  ^ 

where  ac=— 28.493  mg  m-3,  ai=0.014  mg  m  year-1, 
Z?o=55.051  mg  m-3,  bi=— 0.0275  mg  m  3  year-1. 

4.2.  Optimal  spectral  decomposition  method 

The  OSD  method  was  described  in  detail  for  closed 
(Eremeev  et  al.,  1994,  1995;  Ivanov  and  Margolina, 
1996;  Danilov  et  al.,  1996;  Chu,  1999;  Ivanov  et  al., 
2001,  2004;  Chu  et  al.,  2003a, b)  and  semi-closed 
basins  (Chu  et  al.,  2003a, b;  Ivanov  et  al.,  2004).  Here,  a 
brief  description  is  given  to  illustrate  the  special 
features  of  the  OSD  method  for  analyzing  the 
chlorophyll-^  in-situ  measurements  in  the  Black  Sea. 

4.2.1.  Spectral  decomposition 

The  de-trended  Black  Sea  chlorophyll-a  con¬ 
centration  c(x,  z,  t )  at  depth  zk  is  decomposed 


where  V/*  =  3  /3x  +3  /3 y  ,  and  n  is  the  unit  vector 
normal  to  T(z).  The  basis  functions  are 

independent  of  the  data  and  therefore  available  prior 
to  the  data  analysis. 

The  OSD  method  has  two  important  procedures: 
optimal  mode  truncation  and  determination  of  spectral 
coefficients  {Amj.  After  the  two  procedures,  the 
generalized  Fourier  spectrum  (3)  is  used  to  provide 
data  at  regular  grids  in  space  and  time. 

4.2.2.  Optimal  mode  truncation 

The  optimal  mode  truncation  number  (Mopt)  is 
defined  as  the  critical  mode  number  with  the  set  of 
spectral  coefficients  {Am}  least  sensitive  to  observatio¬ 
nal  data  sampling  and  noise.  For  sample  size  of  P  and 
mode  truncation  of  M,  the  spectral  coefficients  {Am} 
are  estimated  by  the  least  square  difference  between 
observed  and  calculated  values  (Menke,  1984), 

d smp  =  J  (A  i ,  .  .  .  ,  A  yf ,  P ,  M  ) 

1  P  (  .  M  \2 

=  pE  cW  “  'S23m(z,t)V^(x,z)  “ >min, 

J=  1  V  m=\  ) 

(5) 

where  the  symbol  denotes  the  estimated  values  at 
location  (x,  z).  For  homogeneously  sampled  data  with 
low  noise  and  without  systematic  error,  the  empirical 
cost  function  Jemp  should  tend  to  0  monotonically  as  M 
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increases  to  infinity.  The  set  of  the  spectral  coefficients 
{Am}  depends  on  the  mode  truncation  M.  Optimal 
estimation  of  {Am}  is  equivalent  to  the  determination 
of Mopt  (Ivanov  et  al,  2001, 2004;  Chu  et  al,  2003a, b). 

A  modified  cost  function  (Vapnik  and  Chervonen- 
kis,  1974a,b), 


demp 

■-\ 

M\ 

K*0 

|  -  ln(l  -  t) 

P 

is  used  to  determine  the  optimal  mode  truncation 
Here,  t  is  the  probability  of  \J— ./empl- '0  as  M 
increases.  It  is  reasonable  to  take  t=0.95  (Mikhal’sky, 
1987).  Previous  studies  on  the  Black  Sea  chlorophyll- 
a  and  radionuclides  Cs-137,  Cs-134  and  Sr-90 
(Eremeev  et  al.,  1994,  1995;  Ivanov  et  al.,  2001) 
show  that  Mopt  is  30-50  for  the  basin-scale  (-300  km) 
variability  and  150  for  the  mesoscale  (-20  km) 
variability.  For  sparse  and  noisy  data,  it  is  difficult 
to  get  reliable  and  stable  estimates  of  all  the  necessary 
spectral  coefficients,  but  the  first  six  spectral  coef¬ 
ficients  A0(zk,t),  Ai(zk,t),. . .,  A5(zk,t )  are  reliable  and 
stable.  This  is  consistent  with  an  earlier  study  (Ivanov 
et  al.,  2001).  Since  the  first  six  modes  represent 
neither  basin-scale  variability  (30-50  modes  required) 
nor  mesoscale  variability  (150  modes  required),  the 
only  physically  realistic  estimation  is  the  horizontally 
averaged  concentration  c(z,t ),  which  is  A0(zk,t). 

4.2.3.  Regularization  procedure 

Determination  of  the  spectral  coefficients  is 
achieved  by  solving  a  set  of  linear  algebraic  equations 
of  {Am(z,t) }  that  are  obtained  from  the  optimization 
procedures  (3)  and  (6).  Due  to  the  high  level  of  noise 
contained  in  the  observations,  the  set  of  algebraic 
equations  is  ill-posed  and  needs  to  be  solved  by  a 
regularization  method  that  requires:  (a)  stability 
(robustness)  even  for  data  with  high  noise,  and  (b) 
the  ability  to  filter  out  errors  with  a-priori  unknown 
statistics.  Numerical  tests  (Ivanov  et  al.,  2001)  on  the 
reconstruction  of  the  Black  Sea  June  climatological 
chlorophyll-a  concentration  demonstrate  the  power¬ 
fulness  of  the  technique  on  the  basis  of  Eq.  (6)  in 
conjunction  with  Tikhonov’s-like  regularization  pro¬ 
cedures  (Tikhonov  et  al.,  1990)  and  the  rotation 
matrix  method  (Chu  et  al.,  2004). 


4.3.  Filling  of  temporal  gaps 

Because  of  the  data  sparseness  and  irregularity 
(Fig.  4),  the  horizontally  averaged  concentration 
c{zk4i)  [i.e.,  AvizkJ)']  is  directly  calculated  only  for 
the  following  nine  months:  January,  March,  April, 
June,  July,  August,  September,  October  (only  for 
the  sea  surface)  and  November.  The  Fourier  series 
expansion  is  used  to  find  the  values  for  other 
months, 

L 

c{zk,ti)  =  a0(zk)  +  ^2  [ai(zk)cos(la>ti) 
i= l 

+  bi(zk)sm(la>ti)\,  (7) 

with  L= 3,  cd=2t[/T,  and  T=  12  months.  The  seven 
spectral  coefficients  a0 ,  aj  and  bj  ( 1=1 ,  2,  3)  are 
calculated  from  the  data  [A0(zk,t)\  using  the  least 
square  principle  with  a  natural  a-priori  constraint, 

c{zk,ti)>  0. 

4.4.  Reconstruction  skill 

Clearly,  reconstruction  of  the  Black  Sea  chloro- 
phyll-a  concentration  from  noisy  and  sparse  in-situ 
data  is  an  ill-posed  problem.  The  OSD  method  has  the 
following  merits. 

(i)  Robustness .  The  reconstructed  monthly  mean 
Black  Sea  chlorophyll-a  concentration  is  robust 
to  variation  of  the  observational  sampling  and 
small  perturbations  of  the  measurements  (pro¬ 
perty  of  any  regularized  estimates  (Tikhonov  et 
al.,  1990).  In  many  cases,  the  regularized 
estimations  are  insensitive  to  large  observation 
errors  (Chu  et  al.,  2004).  For  observations  with 
white  Gaussian  noise,  the  vertical  structure  of 
the  monthly  mean  Black  Sea  chlorophyll-a 
concentration  can  be  reconstructed  even  with 
high  noise-to-signal  ratios  such  as  2-4. 

(ii)  Positiveness .  The  chlorophyll-^  concentration 
must  be  positive.  This  is  a  very  strong 
a-priori  constraint  to  filter  numerous  unreal¬ 
istic  solutions. 

(iii)  Statistical  completeness .  The  spectral  coeffi¬ 
cients  Am(z,t)  tend  to  the  true  values  without 
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any  discontinuity  as  the  observations  increase  to 
infinity  (Vapnik  and  Chervonenkis,  1974a,b). 
This  property  is  automatically  preserved  using 
the  cost  function  (6)  if  no  systematic  error  exists 
in  the  observations. 


spectrum  associated  with  the  measurement  error  is 
represented  by  (Sabel’feld,  1991) 


4.5.  Sources  of  reconstruction  errors 

Three  sources  of  uncertainty  exist  in  reconstructing 
the  chlorophyll-a  concentration:  (a)  measurement 
errors,  (b)  unresolved-scale  errors  caused  by  sparse 
and  irregular  sampling. 

4.5.1.  Measurement  errors 

In  the  in- situ  chlorophyll-a  concentration  data, 
approximately  30%  of  the  measurements  use  the 
standard  SP  method  (Jeffrey  and  Humphrey,  1975), 
and  about  70%  of  the  measurements  use  the  FL 
method  (Lofrus  and  Carpenter,  1971).  After  compar¬ 
ing  the  accuracy  of  the  two  methods  from  R/V 
“Professor  Vodyanitsky”  during  March- April  1999, 
Yunev  et  al.  (2002)  found  that  the  maximum  relative 
difference  of  the  chlorophyll-a  concentrations  meas¬ 
ured  by  the  two  methods  is  less  than  19%.  This  allows 
us  to  collect  all  the  measurements  in  a  single  data  set 
without  any  preliminary  correction. 

Individual  chlorophyll-a  concentration  measure¬ 
ment  error  varies  in  space  (from  station  to  station)  and 
time.  The  chlorophyll-a  signals  are  mixed  with  the 
accessory  pigments  (chlorophyll-/?  and  chlorophyll-c) 
and  the  degradation  products  (pheophytin-a,  pheo- 
phytin -b  and  pheophytin-c)  (Hallegraeff  and  Jeffrey, 
1985;  Moberg  et  al.,  2001).  These  pigments  and 
pheopigments  are  interfering  compounds  in  the  SP 
and  FL  methods  because  they  absorb  and  fluoresce  at 
the  same  wavelengths  as  the  chlorophyll-a  (Lorenzen 
and  Jeffrey,  1980).  However,  the  error  estimation  for 
this  case  is  poorly  documented. 

Measurements  in  the  Black  Sea  (Vedernikov  and 
Demidov,  1997)  and  the  Indian,  Pacific  and  Atlantic 
Oceans  (Trees  et  al.,  1985)  show  that  the  SP  method 
overestimates  the  chlorophyll-a  concentration  by 
about  60%;  and  the  FL  method  underestimates  the 
chlorophyll-a  concentration  by  about  30%.  The  mean 
measurement  error  ranges  between  about  40%  and 
50%  from  station  to  station.  As  a  working  hypothesis, 
the  measurement  error  is  assumed  to  be  Gaussian 
white  noise  (0,  a)  in  space  and  time.  The  noise 


where  M0  (~300)  is  the  spectral  mode  analogue  to  the 
Nyquist  frequency  of  a  regular  discrete  time  series. 
For  two-dimensional  sampling,  the  wavelength 
Lo=MflL  has  the  order  of  “twice  the  averaged  station 
separation”.  Here,  L  is  the  integral  horizontal  scale  of 
the  Black  Sea. 

The  variance  of  the  reconstructed  chlorophyll-a 
concentration  field  is  represented  by  the  Euclidean 
norm  ||^^||OT3-  The  variance  of  noise  in  the  input 
signal  is  estimated  by 


=  PMopt/Mo ~  1 .2  x  If)"2,  (9) 

which  shows  that  the  optimal  mode  truncation  Mopt 
reduces  drastically  the  impact  of  measurement  noise 
on  the  reconstruction  accuracy. 


4.5.2.  Unresolved-scale  errors 

The  second  error  source  is  unresolved  synoptic  and 
mesoscale  motions.  The  Black  Sea  has  strong 
mesoscale  and  synoptic  variability,  such  as  in  the 
Rim  current  on  the  shelves  (Ozsoy  and  Unluata,  1997; 
Oguz  et  al.,  2002a),  that  greatly  affects  the  spatial  and 
temporal  variability  of  chlorophyll-a  concentration. 
Based  on  numerical  simulation  and  analysis  of  the 
SeaWiFS  images,  Oguz  and  Salihoglu  (2000)  show 
that  the  mesoscale  features  can  control  patchiness  and 
intensity  of  chlorophyll-a  concentration  ranging 
between  1  and  3  mg  m-3  during  the  spring  bloom. 
As  the  first  guess,  we  may  parameterize  this  varia¬ 
bility  as  red  noise  with  a  correlation  radius  Rcor 
(Thiebaux  and  Pedder,  1987). 

Twin  experiments  show  weak  sensitivity  of  the 
reconstruction  skill  of  a  scalar  field  in  the  Black  Sea  to 
the  correlation  radius  of  the  unresolved-scale  pro¬ 
cesses  represented  by  the  Gaussian  red  noise  if 
Rcor<2L0  (Ivanov  et  al.,  2001).  The  reconstruction 
error  grows  approximately  linearly  for  large  correla¬ 
tion  scales.  If  the  Rossby  deformation  radius  Rd 
(around  35-40  km  for  the  Black  Sea)  is  used  to 
represent  the  correlation  radius  ( Rcor~Rd ),  then  Rcor  is 
much  smaller  than  2 L0.  Therefore,  the  variability  with 
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such  spatial  scales  should  be  effectively  filtered  by  the 
regularization  procedure. 

It  is  more  difficult  to  evaluate  the  impact  of 
unresolved-scale  errors  on  the  reconstruction  skill 
than  of  measurement  errors.  Usually,  the  unresolved- 
scale  error  variance  is  assumed  to  be  30-50  times  the 
measurement  error  variance  o2  (Thiebaux  and  Pedder, 
1987),  which  leads  to  40-60%  error  variance  distort¬ 
ing  the  input  signal. 

The  regularization  guarantees  an  accuracy  around 
20-30%  of  the  norm  ||^^||^=i  utilizing  the  obser¬ 
vational  data  with  40-60%  error  variance  (Ivanov  et 
al.,  2001;  Chu  et  al.,  2003a, b).  Thus,  20-30%  can  be 
taken  as  the  low  accuracy  limit  in  estimating  A(zk,t ) 
[i.e.,  c(zk,ti)\. 

4.5.3.  Irregular  sampling  errors 

Irregular  sampling  is  also  a  major  source  of 
uncertainty.  The  OSD  method  is  sufficient  to  filter 
out  errors  if  the  characteristic  spatial  scale  of  the 


reconstructing  field  L0  is  twice  the  characteristic  scale 
of  the  observational  gaps  Zgap.  In  estimating  A(zk,t), 
the  reconstruction  procedure  was  proved  efficient 
(Ivanov  et  al.,  2001)  through  reconstructing  horizontal 
mean  Black  Sea  scalar  fields  with  various  sampling 
strategies. 

5.  Seasonal  variability 

Horizontally  averaged  chlorophyll-a  concentra¬ 
tions  in  the  Black  Sea,  c(zk,ti),  are  calculated  using 
the  OSD  method  [i.e.,  A(zk,t ),  Fig.  5a]  and  the  simple 
arithmetic  mean  (in  the  horizontal  and  month)  of  the 
in-situ  chlorophyll-a  concentration  data  (Vedernikov 
and  Demidov,  1997,  hereafter  called  the  VD97 
profiles,  Fig.  5b)  at  11  depths  zk  [k=l,  K ,  ^=11] 
between  the  sea  surface  and  50  m  depth  with  5-m 
interval.  There  are  11  VD97  profiles  for  all  months 
except  December.  The  distinction  between  the  two  is 
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Fig.  5.  Time-depth  cross-sections  of  c(z,t )  from  (a)  the  OSD  reconstructed  field  with  two  blooms  in  February/March  (winter/spring  bloom)  and 
September/October  (fall  bloom),  (b)  the  VD97  profiles  with  only  one  bloom  in  February/March  (winter/spring  bloom). 
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the  two-bloom  structure  (winter/spring  and  fall)  in  the 
reconstructed  field  (OSD  method)  versus  one-bloom 
structure  (winter/spring)  in  the  VD97  profiles.  The 
seasonal  variability  identified  using  the  OSD  method 
is  described  below. 

Sufficient  nutrient  supply,  light  availability,  and 
biological  interaction  might  trigger  a  bloom.  Sverdrup 
(1953)  defined  a  critical  depth  above  which  the  depth- 
integrated  daily  gross  primary  production  equals 
phytoplankton  loss  through  grazing  and  respiration. 
The  critical  depth  is  associated  with  the  light 
availability  and  deepens  in  spring  and  summer. 

During  the  winter  season,  intensive  mixing  pro¬ 
vides  the  upper  layer  with  nutrients  but  the  phyto¬ 
plankton  production  remains  low  because  of  light 
limitation.  The  winter/spring  bloom  occurs  in  Febru¬ 
ary/March  with  a  maximum  value  of  2.25  mg  m"3.  This 
period  is  characterized  by  near  uniform  chlorophyll-a 
distribution  within  the  mixed  layer  (30  m  depth), 
which  is  shallower  than  the  critical  layer  depth. 
Phytoplankton  transport  to  below  the  critical  depth  is 
prevented  and  the  majority  of  the  phytoplankton 
remains  within  the  euphotic  zone.  This  mechanism  in 
combination  with  the  nutrient  abundance  triggers  the 
winter/spring  bloom  which  continues  until  nutrient 
depletion  occurs  in  the  upper  layer. 

At  the  end  of  April,  spring  bloom  terminates  as  a 
result  of  depletion  of  the  nitrate  stock  and/or  increasing 
role  of  the  grazing  pressure  imposed  by  the  growing 
zooplankton  community.  The  near  surface  chloro¬ 
phyll-a  concentration  decreases,  but  the  subsurface 
chlorophyll-a  maximum  exists  almost  all  the  time, 
except  in  November,  with  varying  depth  and  intensity. 
Depth  of  the  chlorophyll-a  maximum  (exceeding  0.6 
mg  m"3)  is  localized  between  15  and  20  m  in  late 
spring/early  summer.  In  early  fall,  it  dips  to  25-35  m 
depth  and  intensifies  by  nearly  a  factor  of  two.  During 
all  of  the  warm  period  from  May  till  September,  the 
phytoplankton  production  in  the  surface  layer  remains 
low  because  of  nutrient  limitation  even  though  the 
light  conditions  are  favorable  for  photosynthesis. 

Towards  the  end  of  September,  the  surface  chlor¬ 
ophyll  concentration  begins  to  increase  again  and 
reaches  its  second  maximum  in  October.  The  fall 
bloom  has  half  the  intensity  of  the  winter/spring  bloom 
and  is  restricted  within  the  upper  10  m  only.  The 
surface  (above  10  m)  and  subsurface  (below  20  m) 
chlorophyll-a  concentration  maxima  (>1 .0  mg  m'3)  are 


isolated  by  relatively  low  concentrations  ranging 
between  0.5  and  1  mg  m"3. 

Like  the  winter/spring  bloom,  the  fall  bloom  is 
mainly  driven  by  the  mixing  processes  in  the  Black 
Sea.  As  vertical  mixing  intensifies  in  the  early  fall,  the 
seasonal  thermocline  starts  to  erode  and  favors  the 
vertical  transport  of  nutrients  from  deep  to  surface 
waters  and  of  phytoplankton  to  depth  with  favorable 
insolation  conditions.  Oguz  et  al,  (2000,  2003)  pointed 
out  that  the  fall  bloom  is  governed  not  only  by  physical 
processes  but  also  by  biological  interactions  within  the 
Black  Sea  environment,  such  as  fall  rebounds  in 
gelatinous  carnivore  populations  and  subsequent 
decreases  in  herbivorous  mesozooplankton  stocks. 

The  two  maxima  (winter/spring  and  fall)  and  one 
minimum  in  the  summer  are  generally  recognized  as  a 
natural  pattern  of  the  phytoplankton  annual  cycle  in 
middle  latitudes  (Longhurst,  1995;  Sathyendranath  et 
al.,  1995).  Physically  driven  (wind-induced,  fronts, 
upwelling)  water  mixing  drives  seasonal  nutrient 
supply  to  the  euphotic  zone,  and  in  turn  causes  the 
primary  production  (Longhurst,  1995).  In  addition  to 
this  “bottom-up”  forcing,  “top-down”  control  by  the 
gelatinous  predators  plays  a  remarkable  role  in  the 
Black  Sea  phytoplankton  distribution  and  evolution 
on  the  basin  scale.  This  was  confirmed  in  numerical 
simulations  using  a  one-dimensional,  vertically 
resolved,  coupled  physical  biochemical  model  (Oguz 
et  al.,  2000). 


6.  Subsurface  chlorophyll-a  maximum 

The  reconstructed  data  (OSD  method)  show  that  the 
subsurface  chlorophyll-a  maximum  is  located  at  40  m 
depth  in  April  and  35  m  depth  in  August  (Fig.  6a), 
occurring  just  after  the  winter/spring  bloom  and  shortly 
before  the  fall  bloom.  The  shoaling  of  the  chlorophyll- 
a  maximum  occurs  in  June  (15  m)  and  in  December- 
February  (15  m)  between  the  two  deepening  events  (in 
April  and  August)  and  is  accompanied  by  a  continuous 
reduction  of  subsurface  maximum  chlorophyll-a  con¬ 
centration  (Fig.  6b).  The  subsurface  maximum  peaks  at 
2.25  mg  m-3  in  February  and  peaks  again  at  1.16  mg 
m-3  in  September-October.  Two  minima  occur  0.40 
mg  m-3  in  July  and  0.58  mg  m-3  in  December. 

The  VD97  profiles  show  that  the  subsurface 
chlorophyll-a  maximum  is  located  at  5  m  depth  in 
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Vedernikov  &  Demidov,  1997 

Fig.  6.  Monthly  mean  (a)  subsurface  chlorophyll-a  maximum  depth,  and  (b)  subsurface  chlorophyll-a  maximum  concentration,  estimated  using 
OSD  (dotted)  and  simple  arithmetic  mean  (Vedernikov  and  Demidov,  1997)  (circle). 


the  winter  (November- January)  and  45  m  in  June  (Fig. 
6a).  The  subsurface  chlorophyll-a  concentration  in  this 
data  set  has  one  maximum  of  2.75  mg  m-3  in  February 
and  one  minimum  of  0.25  mg  m-3  in  May  (Fig.  6b). 

Two  major  discrepancies  are  found  between  the 
reconstructed  and  VD97  profiles:  (a)  the  subsurface 
chlorophyll-a  concentration  in  September-0 ctober  is 
1.16  mg  m-3  in  the  reconstructed  profiles  and  0.50- 
0.55  mg  m-3  in  the  VD97  profiles,  and  (b)  the  depth  of 
the  subsurface  chlorophyll-a  maximum  in  June  is  1 5  m 
in  the  reconstructed  profiles  and  45  m  in  the  VD97 
profiles. 

In  the  fall  (September-October),  the  depth  of  the 
subsurface  chlorophyll-a  maximum  is  25  m  in  the 
reconstructed  profiles  and  30  m  in  the  VD97  profiles, 
both  occur  in  the  upper  mixed  layer  with  favorable 
photic  conditions.  These  conditions  may  favor  the 
existence  of  the  second  maximum  in  the  subsurface 
chlorophyll-a  concentration. 

In  June,  whether  the  true  depth  of  the  subsurface 
chlorophyll-a  maximum  is  shallow  (15  m)  or  deep  (45 
m)  can  be  evaluated  using  the  existing  analysis  and 
data.  Based  on  the  data  collected  during  the  last  20 


years,  Yunev  et  al.  (1999)  show  the  existence  of  four 
major  types  of  chlorophyll-a  concentration  profiles. 
According  to  their  classification,  one  well  advanced 
and  symmetric  maximum  or  two  and  more  maxima, 
which  are  comparable  or  differing  in  magnitude,  are 
characteristic  features  of  the  vertical  distribution  of 
chlorophyll-a  during  summer  and  autumn  when  the 
temperature  stratification  is  a  barrier  to  penetration  of 
additional  nutrients  from  the  bottom  layers  into 
surface  layer.  For  example,  two  in-situ  measured 
chlorophyll-^  concentration  profiles  (Fig.  7)  clearly 
indicate  the  existence  of  a  chlorophyll-a  maximum  at 
nearly  20  m  depth  (Oguz  et  al.,  1999). 

The  location  and  intensity  of  the  subsurface 
chlorophyll-a  maximum  in  the  Black  Sea  are  gov¬ 
erned  by  nutrient  availability  in  the  euphotic  zone 
(Yunev  et  al.,  1999).  The  winter/spring  high  primary 
production  causes  subsequent  depletion  of  the 
nutrients  in  the  surface  layer,  but  not  in  the 
thermocline.  This  increases  the  chlorophyll-a  concen¬ 
tration  in  the  algal  cells  and  population  of  phyto¬ 
plankton  in  the  thermocline,  and  establishes  the 
subsurface  chlorophyll-^  maximum  in  the  stratified 
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C,  mg/m3 


Fig.  7.  Two  observed  summer  vertical  profiles  of  chlorophyll-* 

ocean  (Kirk,  1983).  Thus,  the  subsurface  chlorophyll- 
a  maximum  should  be  located  just  below  the  seasonal 
thermocline.  On  the  other  hand,  the  layer  of  the 
chlorophyll-a  abundance  is  restricted  by  the  upper 
boundary  of  the  permanent  pycnocline.  In  June,  the 
upper  boundary  of  seasonal  thermocline  (permanent 
pycnocline)  is  10  m  (49-52  m)  for  the  central  western 
(Fig.  8a)  and  eastern  (Fig.  8b)  cyclonic  gyres.  The 
estimated  depth  of  the  subsurface  chlorophyll-a 
maximum  in  June  in  the  reconstructed  profiles  (15 
m)  is  more  realistic  than  in  the  VD97  profiles  (45  m). 


7.  Total  chlorophyll-a  amount 

Vertical  integration  of  A(zk,t)  (i.e.,  estimation  of 
the  horizontally  averaged  chlorophyll-a  concentra¬ 
tion)  from  the  sea  surface  to  50  m  depth,  (9 50ft)  is 
taken  as  the  monthly  mean  total  chlorophyll-a  amount 
with  the  annual  mean 

(Cso)  =2  fT Qso(t)dt~ 0.40  g  Chi  m-2,  (11) 

1  Jo 


C,  mg/m3 


in  the  Black  Sea:  (a)  June  20,  1986,  and  (b)  May  20,  1984. 
and  variance 

(ZQ25o)=yj\Q50-(Q5o))2dt~0.01  (g  Chl)2  m-4. 

(12) 

The  total  chlorophyll-a  amount  from  the 

reconstructed  data  (Fig.  9)  has  two  maxima 
occurring  in  February  [1.00  g  (Chl)  m-2]  and 
September/October  [0.48  g  (Chl)  m-2]  and  two 
minima  occurring  in  July-August  and  November- 
December  [0.1  g  (Chl)  m-2].  The  variability  in¬ 
dices  (i.e.,  ratio  of  maximum  to  minimum)  are  9.9 
and  4.9  for  the  two  cycles,  respectively.  The 
monthly  mean  phosphate  and  depth  of  the  photo¬ 
synthetic  layer  are  also  plotted  in  Fig.  9  to  help 
understand  the  phytoplankton  dynamics  using  the 
abiotic  factors. 

The  phosphate  concentration  has  a  maximum  [0.77 
pmol  (P)  l~l]  in  February.  The  available  nutrients 
have  been  transported  to  the  upper  layer  in  the  winter 
due  to  the  intensive  mixing  processes.  With  the 
shoaling  of  the  permanent  pycnocline  to  about  30 
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-O-  upper  boundary  of  the  permanent  pycnocline 
-O-  upper  boundary  of  the  seasonal  thermocline 

Fig.  8.  Upper  boundaries  of  the  permanent  pycnocline  and  seasonal  thermocline  in  the  (a)  western  and  (b)  eastern  cyclonic  gyres. 


m,  the  phytoplankton  remain  in  the  euphotic  zone 
with  the  shallowest  depth  of  30  m.  This  causes  the 
major  bloom  in  February. 


After  the  winter/spring  bloom,  the  upper  layer  total 
chlorophyll-a  amount  reduces  rapidly  to  0.27  g  (Chi) 
m-2  in  May  because  of  nutrient  limitation  despite  a 
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Fig.  9.  OSD  reconstructed  monthly  mean  Q50  with  error  bars. 
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favorable  combination  of  increased  daily  insolation 
and  onset  of  thermal  stratification.  From  May  to  June 
(early  summer),  the  rate  of  this  reduction  slows  down 
possibly  due  to  the  intensification  of  the  subsurface 
chlorophyll-a  maximum,  which  was  discussed  in  the 
previous  section. 

In  July/August,  the  upper  layer  total  chlorophyll-a 
amount  (?5o(h)  starts  to  increase  and  reaches  its 
second  maximum  value  [0.48  g  (Chi)  m-2]  by 
October,  associated  with  the  fall  bloom.  The  strength 
of  the  fall  bloom  is  nearly  half  that  of  the  winter/ 
spring  bloom.  During  that  period,  the  seasonal 
thermocline  starts  to  erode  because  of  the  strong 
surface  atmospheric  forcing,  which  is  followed  by  the 
nutrient  upwelling  and  phytoplankton  rising  to  depths 
of  favorable  insolation  conditions. 


8.  Significance  of  the  fall  bloom 

The  Black  Sea  has  poor  data  coverage  of  the 
chlorophyll-^  concentration  in  February  and  October, 
except  at  the  surface  (Table  1  and  Fig.  4).  But  the 
reconstructed  chlorophyll-a  field  shows  the  blooms 
during  these  months.  The  poor  data  average  may  raise 


reasonable  doubts  about  the  reliability  of  the  recon¬ 
structed  field.  To  show  the  effectiveness  and  stability 
of  this  temporal  interpolation  scheme,  the  seven 
spectral  coefficients  {d0,ahbi\l=l,2,3}  in  Eq.  (7)  are 
calculated  in  the  upper  5-m  layer,  A0(z>—5  m,^), 
from  (a)  all  the  available  data,  and  from  excluding  the 
data  measured  in  (b)  April,  and  (c)  August  with  the 
constraint,  c(zk,ti)> 0  (Fig.  10). 

The  two  blooms  (winter/spring  and  fall)  exist  in  all 
three  reconstructed  sets  and  little  difference  is  found 
among  them.  The  difference  is  very  small  comparing 
(a)  and  (b),  and  is  less  than  12%  for  the  winter/spring 
bloom  and  less  than  30%  for  the  fall  bloom  between 
using  (a)  and  (c).  The  ratio  of  the  first-to-second 
maximum  chlorophyll-a  concentration  is  1.8,  2.2  and 
2.3  for  the  three  data  sets.  The  stability  of  the 
reconstructed  chlorophyll-a  seasonal  variability 
shows  the  power  of  the  Fourier  series  expansion  to 
effectively  obtain  the  seasonal  variability  (two 
blooms)  using  temporally  interrupted  data.  Although 
it  has  been  identified  using  the  SeaWiFS  chlorophyll 
data  (Oguz  et  al.,  2002a, b),  the  fall  bloom  (Septem- 
ber-October)  in  the  Black  Sea  is  first  detected  in  this 
study  using  the  in-situ  measurements  with  insufficient 
sampling  in  this  season. 


Fig.  10.  Sensitivity  of  OSD  reconstruction  of  the  upper  layer  (within  5  m  depth)  chlorophyll-a  concentration  to  temporal  interruption  using  all 
available  data  (solid  curve),  data  excluding  April  (dashed  curve),  and  data  excluding  August  (dotted  curve). 
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9.  Sensitivity  to  sampling  size 

The  data  for  VD97  profiles  were  collected  at  16 
expeditions  from  1978  to  1992.  The  data  for  Yunev  et 
al.’s  (2002)  estimation  contain  about  1000  observa¬ 
tions  in  the  deep  region  of  the  Black  Sea  in  1964, 
1973,  1978,  and  from  1980  to  1996.  The  data  for  this 
study  contain  about  1100  stations  from  26  cruises  in 
1980-1995  (around  70%  of  all  the  existing  chlor- 
ophyll-a  concentration  data)  and  coincide  partially 
with  the  two  earlier  studies. 

Since  not  all  the  Black  Sea  chlorophyll-a  obser¬ 
vational  data  are  available  for  this  study  (as  well  as  for 
others),  sensitivity  of  the  reconstruction  to  the 
sampling  size  of  our  data  set  was  investigated.  The 
results  show  that  a  reduction  of  10-15%  in  sampling 
size  affects  the  estimation  of  the  monthly  mean 
chlorophyll-a  concentration  slightly.  However,  the 
sampling  size  cannot  be  reduced  any  more  than  10- 
15%,  since  further  reduction  will  impair  the  filtering 
capabilities. 

10.  Effect  of  shelf  data 

Vedernikov  and  Demidov  (1997)  and  Yunev  et  al. 
(2002)  utilized  the  in-situ  data  for  the  deep  basin  of 
the  Black  Sea  only  (with  water  depth  larger  than  200 
m).  In  the  present  study,  we  use  the  in-situ  data  on  the 
northwestern  shelf  of  the  Black  Sea  along  with  the 
deep-basin  data.  Since  the  purpose  of  this  research  is 
to  estimate  the  seasonal  variability  of  horizontal  mean 
chlorophyll-a  concentration,  it  is  necessary  to  assess 
the  effect  of  the  shelf  data. 

The  effect  of  the  shelf  data  on  the  seasonal 
variability  of  the  horizontal  mean  chlorophyll-^ 
concentration  is  significant  using  the  simple  arith¬ 
metic  mean  (the  VD97  approach),  and  insignificant 
using  the  OSD  method.  This  is  because  the  truncated 
expansion  (Eq.  (3))  contains  various  horizontal  scales 
of  the  chlorophyll-a  field  with  the  mode-0  represent¬ 
ing  the  horizontal  mean  field.  The  first  30-50  modes 
represent  the  basin-scale  (300  km)  variability,  and  the 
first  150  modes  represent  the  mesoscale  (20  km) 
variability  including  the  variability  on  the  northwest¬ 
ern  Black  Sea  shelf.  Low-order  truncation  (at  mode- 5 
in  this  study)  automatically  removes  the  effects  of 
smaller  scale  (smaller  than  the  basin-scale)  processes 


including  the  mesoscale  processes  on  the  northwest¬ 
ern  Black  Sea  shelf. 

Let  Qd+s  and  Qd  be  the  horizontally  averaged 
vertical  integrations  (0-5  m  depth)  of  the  chlor¬ 
ophyll-#  concentration  over  the  whole  domain  and 
the  deep-basin  only.  The  ratio  QJQd+5  represents 
the  effect  of  the  shelf  data  with  little  effect  if  Qdl 
Qd+ 5~1,  and  huge  effect  if  QJQd+ 5~0.  The  ratio  for 
the  September  chlorophyll-#  field  is  0.89  using  the 
OSD  method  and  0.30  using  the  simple  arithmetic 
method.  This  shows  the  capability  of  the  OSD 
method  to  filter  processes  with  scales  smaller  than 
the  basin-scale. 


11.  Conclusions 

(1)  Noisy,  sparse,  and  irregularly  distributed  chlor¬ 
ophyll-#  concentration  data  collected  in  the 
Black  Sea  (including  shelf  waters)  between 
1980  and  1995  are  analyzed  using  the  OSD 
method.  The  reconstruction  results  present  the 
depth-resolved  annual  variations  of  the  hori¬ 
zontally  averaged  chlorophyll-#  concentration 
within  the  upper  50-m  layer. 

(2)  The  reconstructed  chlorophyll-#  concentration 
shows  the  existence  of  two  blooms  with  a 
strong  winter/spring  bloom  [up  to  2.00  mg  (Chi) 
m-3]  in  Lebruary/March  and  a  less  strong  fall 
bloom  [around  1.00  mg  (Chi)  m-3]  in  Septem¬ 
ber/October.  The  subsurface  chlorophyll-#  con¬ 
centration  also  has  a  bi-modal  temporal 
structure:  two  maxima,  with  the  first  maximum 
of  2.25  mg  m-3  in  Lebruary  and  the  second 
maximum  of  1.16  mg  m  in  September- 
October,  and  two  minima,  with  the  first  mini¬ 
mum  of  0.40  mg  m-3  in  July  and  the  second 
minimum  of  0.58  mg  m-3  in  December. 

(3)  The  depth  of  the  subsurface  chlorophyll-# 
maximum  has  two  maxima  and  two  minima. 
The  first  maximum  depth  is  40  m  in  April, 
occurring  just  after  the  winter/spring  bloom,  and 
the  second  maximum  depth  is  35  m  in  August, 
appearing  shortly  before  the  fall  bloom.  The 
minimum  depth  is  15  m,  occurring  in  Decem- 
ber-Lebruary  and  June.  The  noticeable  rise  of 
the  chlorophyll-#  enriched  zone  in  June 
between  the  two  deepening  events  (in  April 
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and  August)  is  accompanied  by  a  near-contin¬ 
uous  reduction  of  upper  50-m  layer  chlorophyll- 
a  concentration. 

(4)  The  winter/spring  bloom  is  the  most  prominent 
signal  in  the  in-situ  measurements  of  chloro- 
phyll-a  concentration  and  was  identified  using  a 
simple  arithmetic  averaging  procedure.  This 
bloom  is  caused  by  sufficient  nutrients  being 
transported  to  the  upper  layer  during  the  winter 
due  to  intensive  mixing  processes  and  a  mixed 
layer  depth  shallower  than  the  critical  depth 
during  the  winter/spring  period  due  to  a 
decrease  of  the  mixed  layer  depth  and  an 
increase  of  the  critical  depth. 

(5)  The  fall  (September/October)  bloom  was 
detected  from  the  SeaWiFS  images  and  repro¬ 
duced  in  coupled  physical-biological  models, 
but  cannot  be  identified  from  the  in-situ  meas¬ 
urements  using  the  simple  arithmetic  averaging 
procedure.  This  study  using  the  OSD  method  is 
the  first  time  to  identify  the  fall  bloom  in  the 
Black  Sea  from  the  in-situ  data. 

(6)  The  OSD  method  is  a  useful  alternative  for  the 
popular  optimal  interpolation  method  in  marine 
data  analysis  with  unknown  first  guess  field  and 
autocorrelations  and  high  noise-to-signal  ratio. 

(7)  Mechanisms  explaining  the  existence  of  two 
blooms  and  the  characteristics  of  subsurface 
chlorophyll-a  maximum  in  the  Black  Sea 
should  be  further  investigated  utilizing  addi¬ 
tional  information  on  various  abiotic  factors  that 
affect  the  phytoplankton  dynamics. 
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age  over  the  eastern  than  the  western  Black  Sea 
(Fig.  3).  Why  does  the  reconstructed  horizontally 
averaged  chlorophyll-a  field  represent  the  whole 
Black  Sea  when  the  east  Black  Sea  is  poorly 
sampled? 

The  spectral  coefficients  {Am}  in  the  Fourier  series 
(3)  can  only  be  approximately  estimated  due  to  the 
existence  of  noise, 

Am  Am  -f  SAm , 

where  {Am}  are  the  “true”  coefficients,  and  {SAmj  are 
errors.  The  total  error  is  computed  by 


Q=  £(<H*)2, 

m= 1 

which  should  not  be  greater  than  the  tolerance  S2, 
Q<S2. 

If  the  basis  functions  { }  are  the  eigenfunctions  of 
the  Laplacian  operator,  the  Fourier  series  (3)  is 
uniformly  convergent  (Morse  and  Feshbach,  1953) 
and  the  following  theorem  exists. 

Theorem.  If  the  Fourier  series  (3)  converges  uni¬ 
formly  and 

<S2M(<5)^0  as  <5— >0 

then  for  any  spatial  point  x,y^S, 


lim 

c5^0 


M{S) 


c(x,y,z,t)  —  y ^Am(t)Wm(x,y,z) 


m= 0 


-*0, 


(Al) 


which  shows  that  the  truncated  mode  number  M(5) 
always  exists  and  determination  of  M(d)  strongly 
depends  on  the  observation  error.  For  homogeneously 
convergent  Fourier  series,  the  spectral  coefficients 
{Am}  may  also  be  determined  from  a  sub-area, 


* 

So 


/  M  \  2 

I  c  —  Am  ]  ds~ ‘min, 

V  *= i  / 


S0cS, 


(A2) 


Appendix  A.  Representation  of  basin-scale  field 
using  sub-area  data 

The  spatial  distribution  of  the  chlorophyll-a 
observations  is  heterogeneous  with  much  less  cover- 


where  S  is  the  total  area  and  S0  is  an  arbitrary  area 
where  the  data  are  collected.  This  theorem  was  proved 
by  Tikhonov  et  al.  (1990)  and  is  only  valid  if  the  basis 
functions  { Wm  }  are  the  eigenfunctions  of  the  Laplace 
operator.  It  may  not  be  valid  if  the  basis  functions 
{ Wm}  are  empirical  orthogonal  functions  because 
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they  are  not  homogeneously  convergent.  For  accurate 
observations  (without  noise  and  systematic  error),  it  is 
sufficient  to  use  sub-area  data  to  reconstruct  a  field  for 
the  entire  region. 

Two  remarks  should  be  noted  before  using  this 
theory.  First,  observational  data  are  noisy  and 
usually  obtained  from  non-uniform  sampling.  The 
spectral  coefficients  are  approximately  estimated 
and  the  reconstructed  field  may  not  be  accurate. 
The  farther  the  location  from  the  observational  sub- 
area,  the  larger  is  the  reconstruction  error  (large  or 
even  negative  values  of  concentration)  may  occur. 
Second,  the  spatial  scale  of  the  field  L  determines 
the  number  of  modes  M  required  for  mode 
truncation.  Since  M(5)  in  Eq.  (Al)  depends  func¬ 
tionally  on  S,  it  is  true  that  M>M  for  high  noise  and/or 
non-uniform  sampling.  That  is  to  say  that  the  mode 
truncation  may  not  be  totally  determined  from 
observations.  If  the  truncation  is  at  the  mod e-M,  the 
reconstructed  field  has  noticeable  distortion  in  sub- 
areas  with  no  observation. 

Ivanov  et  al.  (2001)  use  Black  Sea  twin  experi¬ 
ments  to  examine  the  effect  of  noise  and  sampling 
strategy  on  the  reconstruction  quality  using  four 
nondimensional  ratios:  (1)  noise-to-signal  ratio  ( r\x ) 
computed  for  areas  containing  observations;  (2)  mode 
number-to-observational  number  ratio  (rj2),  (3)  recon¬ 
struction  error-to-signal  ratio  for  areas  containing 
observations  (E i),  and  (4)  reconstruction  error-to- 
signal  ratio  for  areas  containing  no  observations^). 
Here,  rfi  represents  the  noise,  r\2  the  sampling 
strategy,  and  (E u  E2)  the  reconstruction  errors. 

For  a  Black  Sea  basin-scale  scalar  field,  30-50 
modes  should  be  used.  If  the  measurement  errors 
are  represented  by  Gaussian  white  noise  with 
rj  i<l,  and  rj2<0.l,  the  reconstruction  errors  are 
estimated  as 

£i<0.05,  E2<  0.25. 

For  the  same  conditions,  but  assuming  Gaussian 
red  noise  with  the  correlation  radius  equivalent  to 
the  Rossby  radius  of  deformation  (40  km),  the 
reconstruction  errors  are  estimated  as 

(Eu  E2)  cz.  0.20  —  0.25. 

These  estimations  show  that  the  reconstructed 
horizontally  averaged  chlorophyll-a  field  can  repre¬ 


sent  the  whole  Black  Sea  when  the  east  Black  Sea  is 
poorly  sampled. 
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